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We develop a new method for deconvolving the smearing effect of the survey window 
in the analysis of the galaxy multipole power spectra from a redshift survey. This method 
is based on the deconvolution theorem, and is compatible with the use of the fast Fourier 
transform. It is possible to measure the multipole power spectra deconvolved from the 
window effect efficiently. Applying this method to the luminous red galaxy sample of the 
Sloan Digital Sky Survey data release 7 as well as mock catalogues, we demonstrate how the 
method works properly. Using this deconvolution technique, the amplitude of the multipole 
power spectrum is corrected. Besides, the covariance matrices of the deconvolved power 
spectra get quite close to the diagonal form. This is also advantageous in the study of the 
BAO signature. 

§1. Introduction 

The power spectrum is one of the most fundamental tools to characterize the 
large-scale clustering of galaxies (e.g.,[T])). The role of the power spectrum analysis in 
cosmology is becoming more and more important^ Especially, the power spectrum 
of galaxies plays a key role for extracting the baryon acoustic oscillations (BAO),®^ 
which is essential for the dark energy research. The quadruple power spectrum is 
also useful for measuring the redshift-space distortion, which plays a vital role of 
testing gravity on the scales of cosmology.® - ® 

The method developed by Feldman, Kaiser and Peacock (hereafter FKP, [9])) 
is often adopted to measure the power spectrum.® It is known that the power 
spectrum measured using the FKP method, due to finiteness of the survey volume, 
is convolved with a window function. This convolution is expressed schematically 



as. 



PC ° nV(k) = (2^)3 / d3k ' P ^')W(k - k'), (1-1) 



where k is the wavenumber vector, W(k — k')/(27r) 3 describes the survey window, 
and we refer to P conv (k) as the 'convolved' power spectrum. In the limit of an 
infinitely large survey volume, Wi}s. — k')/(27r) 3 approaches the three dimensional 
delta function, hence P conv (k) becomes the same as the true power spectrum -P(k). 

In a realistic case, however, a survey volume is finite. The window convolution 
modifies a measured power spectrum compared with the true power spectrum. It 
has the following two aspects: (i) change of the amplitude of the power spectrum, 
(ii) mixing of modes with different wavenumbers. These two effects are influential in 
the power spectrum analysis for extracting the BAO and in testing gravity on the 
scales of cosmology. When comparing the convolved power spectrum with theoretical 



typeset using PTPTeX.cIs (Ver.0.9) 



2 



T. Sato, G. Hiitsi, and K. Yamamoto 



predictions, we need to take the window effect into account following Eq. (|l-ip . even 
though it is quite a time consuming process.^ 

If we could directly measure a deconvolved power spectrum from a galaxy data, 
it could be compared with theoretical models without including the window effect 
thus speeding up the analysis significantly. In the present paper, for the first time, 
we develop a scheme to measure the deconvolved multipole power spectra. We apply 
it to the Sloan Digital Sky Survey (SDSS) luminous red galaxy (LRG) sample from 
the data release (DR) 7 as well as to 1000 mock catalogues mimicking the LRG 
sample. This paper is organised as follows: In section 2, we present a new method 
to deconvolve the window function from the FKP estimator, including a brief review 
for the FKP estimator for the power spectrum. In section 3, we apply the method 
to the SDSS LRG DR7 as well as to 1000 mock catalogues, and demonstrate how it 
works in the multipole power spectrum analysis, by comparing the convolved power 
spectrum and the deconvolved power spectrum. Throughout this paper, we use units 
in which the velocity of light equals 1, and adopt the Hubble parameter Hq = WOh 
km/s/Mpc with h = 0.7. 

§2. Formulation 

Here we provide basic formulas for the deconvolution of the window function. 
Taking the Fourier transformation of Eq. (|l-ip , we have 



J al 3 ke- ik s P conv (k) = -^-3 J d 3 k'e- lk ' s P(k') J d 3 ke~ iks W (k) 



• (2-1) 



This can be done exactly as long as the k-space is infinitely large. Then, we have 



/ 



d 3 k'e- ty - s P conv (k') 

d 3 k'e- ik ' s P{k') = (2vr) 3 ^ 7 , (2-2) 

/ d 3 k"e~ lk " s W{k") 



whose inverse transformation yields 

! d 3 k'e- ik '- s P conv (k') 

P(k) = / d 3 se iks J— . (2-3) 

J / d 3 k"e- ik "- s W(k") 

For a discrete density field, we need to subtract the shot-noise contribution. 
Following the FKP method, the estimator of the convolved power spectrum is 

P (k) = [f#snW( S )] (1 + a)5 °' (24) 

where n g (s) is the number density field of galaxies whose mean number density is 
n(s), n rn d(s) is the density field of a random sample whose mean number density is 
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a n(s), ip(s) is the weight factor, s denotes the three dimensional coordinate of the 
redshift-space, and we defined 



J d 3 sn(s)tp 



s 



S ° f d 3 sffl (s)?/; 2 (s)' (2 ' 5) 

which describes the shot-noise contribution. The random catalogue is a set of random 
points without any correlation, which can be constructed through a random process 
by mimicking the selection function of the galaxy catalogue. Note that a is the 
parameter of the random catalogue, for which we adopt the value of 0.01, in the 
present paper. Similarly, we may adopt the estimator for the window function 



/ d 3 san rnd (s)^( 



s e 



;-k|2 



Note that the denominator f d 3 sfi 2 (s)ip 2 (s) cancels out when we substitute Eqs. (f!F 
(|2-5p and (|2-6p into Eq. (|2-3p . Practically, the computation of iy(k) must be done 
carefully, because it might become zero when we choose a box size for the fast Fourier 
transform too large. One can avoid this problem by properly choosing the box size 
in which sample galaxies are distributed. 

The monopole power spectrum is computed from Eq. (j2-3fl . 

1 t j3 



Vk 



P (k) = — / d 6 kP(k), (2-7) 



where Vk is the volume of a shell in the Fourier space, whose mean radius is k. Within 
the small-angle approximation (distant observer approximation), the estimator of the 
higher multipole power spectra might be taken as 



Pe(k) = ^r f d 3 kP(k)£e(e-k), (2-f 



where e denotes the unit vector that points to line of sight direction, £e(fi) is the 
Legendre polynomial, and k = k/|k| is the unit wavenumber vector. Note that the 
line of sight direction is approximated by one direction e in this approach. We refer 
to the multipole power spectrum P^{k) with Eq. (|2-3|) substituted into the right hand 
side of (|2-8p as the 'deconvolved' multipole power spectrum. Note that our definition 
of the multipole spectrum Pi{k) is different from the conventional one by the factor 
21 + i.EJ.03 

In the ref. IT^|) , the FKP method is generalized to measure higher order multipole 
power spectra, which is useful in quantifying the redshift-space distortions (cf. fT2|) , 
[T3])). There are differences between the previous paper-^ 1 and the present one in the 
estimators for the higher multipole spectra. In Table I, the differences of the two 
approaches are summarized. In the previous paperpJ the line of sight direction was 
defined for each pair of galaxies, while in the present paper the line of sight direction 
is approximated by only one direction e. In this sense, the approximation of the 
present paper is limited to a set of galaxies within a narrow survey area. A similar 
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approach was adopted in the references I15 p . [l~6|) . In this approach, because of the 
limitation of the small-angle approximation, a full sample of large survey area must 
be divided into subsamples of narrow survey areas. Because of the division of the 
full sample, the window effect could be very influential, but the advantage of this 
approach is that one can use the fast Fourier transform. 





Present paper 


Reference 11411 


Line of sight direction 


approximated by 


computed for each 




one direction e 


pair of galaxies 


Division of galaxy sample 


necessary 


not necessary 


Window effect 


significant 


not significant 


Use of the fast Fourier transform 


possible 


impossible 



Table I. Comparison of the two approaches 



§3. Application 

In this section we apply our method to the SDSS LRG sample of DR7. Our 
LRG sample is restricted to the redshift range z = 0.16 - 0.47. In order to reduce 
the sidelobes of the survey window we remove some non-contiguous parts of the 
sample, which leads us to 7150 deg 2 sky coverage with the total number N = 100157 
LRGs. The data reduction is the same as that described in the references l3|). [TT|) . 
In our power spectrum analysis we adopt the spatially flat lambda cold dark matter 
model distance-redshift relation with f2 m = 0.28, and set ip = 1. 

We divide the full sample into subsamples, in which galaxies are distributed 
within narrow area. This is necessary to compute the quadrupole power spectrum 
in our method with the fast Fourier transform. In the present paper, we divide the 
full sample into 18 subsamples whose area and numbers of galaxies are almost equal, 
397 square degrees, or into 32 subsamples whose mean area is 223 square degrees. 

3.1. Shape of the multipole power spectrum 

First, we focus on the shapes of the monopole and the quadrupole power spectra. 
Figured] shows the convolved (left panels) and deconvolved (right panels) power spec- 
tra from the SDSS LRG sample. The upper left panel plots the convolved monopole 
spectrum multiplied by the wavenumber, /ci-Q° nv (k). Here the following two cases 
are shown: (i) the full sample divided into 18 subsamples with each piece spanning 
almost 397 square degrees, (ii) the full sample without division into subsamples. 
Due to significantly different effect of the window convolution in these two cases, the 
amplitude of the two spectra are very different. The upper right panel is the same 
as the upper left panel, but is the deconvolved power spectrum kPo(k). Note that, 
in this panel, the amplitudes of the two deconvolved spectra are almost the same. 
The lower panels plot kP2(k). The left panel is the convolved spectrum, while the 
right panel is the deconvolved spectrum. For the quadrupole power spectrum, we 
present the two cases with different divisions: chunks with mean survey area 397 
square degrees, and 223 square degrees, respectively. One can see the same features 
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as those for kPo(k). The results do not depend on the division of the full sample. 
The error-bars are obtained by computing the standard deviations of 1000 mock 
catalogues (A), explained in the next. 

Figure[2]is the same as figure[TJ but shows the average from 1000 mock catalogues 
(A), which are obtained by the following three steps, whose details are described in 
the reference [3]) . First, we generate the density field using an optimized 2nd order 
Lagrangian perturbation calculation. Second, we take the Poisson sampling of the 
generated density field, adjusting the clustering bias and the number density so as 
to agree with the observed LRG sample. Third, we apply the radial and angular 
selection function, and extract the final catalogues. We call these mock catalogues 
(A) . Figure [2] clearly shows that the amplitudes of the two deconvolved spectra are 
almost the same. 




k[h'Mpc-'] 

Fig. 1. Comparison of the convolved fcP f conv (fc) (left panel) and deconvolved power spectra kP l (k) 
(right panel). The upper panels plot the monopole I = 0, while the lower ones the quadrupole 
1 = 2. For the monopole power spectrum, the two cases are presented: One is the case when the 
full sample is divided into 18 subsamples and the other is the case without division of the full 
sample. For the quadrupole spectrum, the two different cases of the division of the full sample 
are presented, in which the mean survey areas are 397 square degrees and 223 square degrees, 
respectively. The error-bars are obtained by computing the standard deviations of 1000 mock 
catalogues (A). 

Figure [3] demonstrates how the deconvolution recovers the original power spec- 
trum properly. To this end, we buid other catalogues (B) and (C). The cosmological 
parameters are the same as those of the mock catalogues (A). To buid the cat- 
alogues (B), the first step is the same as those for the mock catalogues (A), we 
generate the density field using an optimized 2nd order Lagrangian perturbation 
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Fig. 2. Same as figure [T] but using the average from 1000 mock catalogues (A). The error-bars are 
the standard deviations. 



calculation. Then, we take the Poisson sampling of the generated density field, 
adjusting the clustering bias and the number density. Here, we assumed the homo- 
geneous mean number density with n = 8 x 10~ 5 (/i~ 1 Mpc)~ 3 in the square region 
2560/i~ 1 Mpc x 2560/i~ 1 Mpc x 1280/i _1 Mpc. The catalogues (B) are obtained by 
applying the same selection of the survey region as the SDSS LRG sample. Hence, 

(B) has the same survey region as (A) but with the constant mean number density. 
We also consider the catalogues (C), which are the same as (B) but without the 

selection of the survey region. Hence, each catalogue of (C) has the square region 
2560/i~ 1 Mpc x 2560/i~ 1 Mpc x 1280/i _1 Mpc. To obtain the redshift-space catalogues, 
however, for the catalogues (C), we assume an observer at an infinite distance away, 
in order to guarantee the validity of the distant observer approximation. Thus the 
catalogues (C) are ideal samples that assume the large volume and the validity of 
the distant observer approximation, which could be expected to give us the original 
power spectrum. 

Figure [3] plots the convolved spectra (left panels) and the deconvolved spectra 
(right panels) from the catalogues (B). The upper panels are the monopole, while the 
lower panels the quadrupole spectra. In each panel, the convolved spectrum (dotted 
curve) and the deconvolved spectrum (solid curve) from the ideal catalogues (C) 
are shown for comparison. These two curves are almost same, hence the convolved 
spectrum and the deconvolved spectrum are almost same, because the catalogues 

(C) are ideal samples with the large volume. In the right panels of Fig. [3l the good 
agreement of all the spectra indicates that our deconvolution recovers the original 
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Fig. 3. Same as figure [T] but using the average from the 1000 catalogues (B). The solid curve and 
the dotted curve are the deconvolved spectrum and the convolved spectrum from the catalogues 
(C). 



power spectrum properly. However, one can see small deviation in the quadrupole 
spectrum between the result from the catalogues (B) and (C), in the lower right panel. 
As the deviation is larger at small k, one possible reason of this small deviation is 
the limitation of the distant observer approximation for the catalogue (B). 

3.2. Covariance matrix 

In the following we discuss the effect of the deconvolution on the covariance 
matrix, which is defined by 

Cw(h, kj) = (APeik^APeikj)) . (3-1) 

We obtain the covariance matrix Cw using 1000 mock catalogues. The correlation 
matrix, which describes the correlation of the errors between different wavenumbers, 
is defined 

n i \ C w (ki,kj) 

r U /(kj, kj) = (3-2) 

^/Ca{ki,ki)Ci: V {kj,kj) 

Figure 0] shows the correlation matrices of the convolved power spectra for £ = 
£' = (left panel), I = £' = 2 (center panel) and I = 0,f = 2 (right panel), 
respectively, from the mock catalogue (A). Here, the full sample is divided into 
18 subsamples, whose mean survey area is 397 square degrees. On the other hand, 
Figure [S] shows the correlation matrices of the deconvolved power spectra for t = £' = 
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(left panel), £ = £' = 2 (center panel) and £ = 0, £' = 2 (right panel) respectively, 
from the mock catalogue (A). The narrower subsample has a window function with 
a broader width, which makes the power spectrum at different wavenumbers more 
strongly correlated. One can see that the errors of the power spectra at different 
wavenumbers are correlated in figure HI while practically de-correlated in figure [5j 
The off-diagonal components of figure [5] are much smaller than those of figure HI 




Fig. 4. The correlation matrices defined by Eq. (|3-2[) . obtained from 1000 mock catalogs: roo(fe, k') 
(left panel), r22(k,k') (center panel) and ro2(k,k') (right panel), respectively. Here we used the 
convolved power spectra of the subsamples with the mean survey area 397 deg 2 . 




Fig. 5. Same as figure [4] but with the deconvolved power spectra. 

3.3. Baryon acoustic oscillation 

In this subsection, we discuss the effect of deconvolution on the BAO signature. 
In general, the BAO signature in a convolved power spectrum is smoothed compared 
with the original power spectrum. Then, it is expected that the BAO signature of 
the deconvolved power spectrum has a larger amplitude compared to that of the 
convolved power spectrum. 

Figure [6] shows the BAO signature of the convolved and deconvolved power 
spectra from the LRG sample (left panels) and from 1000 catalogues (B) (right 
panels), respectively. In each panel, the solid curve and the dotted curve is the BAO 
signature from the deconvolved spectrum and the convolved spectrum, respectively, 
obtained with 1000 ideal catalogues (C). The dashed curve is the theoretical curve 
in the linear theory. 

The results labelled by (1) and (2) are obtained using the full sample without 
division into subsamples. (1) is from the deconvolved power spectrum Po(k), while 
(2) is the convolved spectrum Pg onv (A:). The BAO signature is obtained using the 
no- wiggle power spectrum multiplied by a function of the form a + bk + ck 2 + d/k 2 , 
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where a, b, c, and d are fitting parameters (see ref]18p for details). The errors in 
the left panels are large, but the following features can be seen in the right panels 
(the results from 1000 mock catalogues (B)). The difference between (1) and (2) is 
small, which means that the window effect is not very influential in this case because 
the survey volume of the full SDSS LRG sample is large. The difference between 
(1) and the solid curve and the dotted curve is also small. This means that the 
BAO signature of the deconvolved spectrum (1) well agrees with the original BAO 
signature. 

The results labelled by (3) and (4) are obtained using 18 subsamples divided 
from the full sample. (3) is from the deconvolved power spectrum, while (4) from the 
convolved one, respectively. The difference between (3) and (4) is not negligible at the 
peaks and the troughs due to the window effect, while the difference between (3) and 
the solid curve and the dotted curve is small. This means that the deconvolution is 
successful and recovers the BAO signature, although the window effect is substantial, 
and degrades the BAO signature when the survey volume is small. Practically, we 
need not divide a full sample into subsamples to obtain the BAO signature from 
the monopole spectrum. Then, the window effect on the BAO analysis seems to be 
small compared with errors. However, in a future precise measurement of the BAO 
signature, the prescription for the window effect might become important. 

§4. Summary and conclusions 

Large galaxy surveys are promising tools for the future dark energy research. In 
a redshift survey, the power spectrum analysis is very fundamental for extracting the 
BAO signature. Similarly, tests of the general relativity on the scales of cosmology are 
becoming important related to this topic in cosmology. This is because a modification 
of the gravity theory is an alternative way to explain the cosmic acceleration, instead 
of introducing the dark energy component. The key to distinguish between the dark 
energy and the modified gravity is the evolution of the cosmological perturbations. 
The redshift-space distortion caused by the peculiar motions of the galaxies, provides 
us a useful chance to test gravity on the scale of cosmology.®'^ 1 Here the multipole 
power spectrum plays a key role. 

In this paper, we developed a new method to measure the deconvolved power 
spectra. By applying it to the SDSS DR7 LRG sample as well as to the mock 
catalogues, we demonstrated that the scheme works well. By the deconvolution, we 
can get the power spectra whose amplitudes are properly recovered. This will be 
essential in measuring the redshift-space distortions. The covariance matrices of the 
deconvolved power spectra are practically diagonal. This is also useful in the study 
of the BAO signature. Our method matches with the fast Fourier transform, which 
saves computation time significantly. 

As mentioned at the end of the subsection 3.1, the limitation of the distant ob- 
server approximation might cause a small deviation of the quadrupole power spec- 
trum when compared with theoretical prediction derived on the basis of the distant 
observer approximation. For a precise treatment, the framework of the spherical 
redshift-space distortion is necessary.^ >E3 However, this problem might be softer 
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Fig. 6. Baryon acoustic oscillation signature of the convolved and deconvolved power spectra from 
the LRG sample (left panels) and 1000 mock catalogues (B) (right panels), respectively. In each 
panel, the solid curve and the dotted curve is from the deconvolved spectrum and the convolved 
spectrum of the ideal catalogues (C). (1) and (2) are obtained using the full sample without 
division into subsamples. (1) denoted by A is from the deconvolved power spectrum Po(k), 
while (2) denoted by + is from the convolved spectrum Po° nv (fc). The results labelled by (3) 
and (4) are obtained using 18 subsamples divided from the full sample. (3) denoted by □ is 
from the deconvolved power spectrum, while (4) denoted by x is from the convolved spectrum, 
respectively. The dashed curve is the theoretical curve in the linear theory. 



for a galaxy sample at higher redshift, because a volume within a small angular area 
can be large. 
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